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ABSTRACT 

We investigate the correlation between the supermassive black holes (SMBHs) mass 
(Afbh) and the stellar velocity dispersion (ct*) in two types of host galaxies: the early-type 
bulges (disk galaxies with classical bulges or elliptical galaxies) and pseudobulges. In the 
form log(Mbh/MQ) = a + /31og(a-*/200 km s"^), the best-fit results for the 39 early-type 
bulges are the slope /3 = 4.06 ± 0.28 and the normalization a = 8.28 ± 0.05; the best-fit 
results for the 9 pseudobulges are /3 = 4.5 ± 1.3, a = 7.50 ± 0.18. Both relations have 
intrinsic scatter in log Afbh of <0.27 dex. The Mbh-c* relation for pseudobulges is different 
from the relation in the early-type bulges over the 3(t significance level. The contrasting re- 
lations indicate the formation and growth histories of SMBHs depend on their host type. The 
discrepancy between the slope of the Mbh-c* relations using different definition of velocity 
dispersion vanishes in our sample, a uniform slope will constrain the coevolution theories of 
the SMBHs and their host galaxies more effectively. We also find the slope for the "core" 
elliptical galaxies at the high mass range of the relation appears steeper (/? ~ 5-6), which may 
be the imprint of their origin of dissipationless mergers. 

Key words: black hole physics - galaxies: bulges - galaxies: formation - galaxies: funda- 
mental parameters - galaxies: nuclei. 



1 INTRODUCTION 

Supermassive black holes (SMBHs) are ubiquitous in the center of 
elliptical galaxies and bulges of disk galaxies. The masses (Afbh) 
of the nearby mostly quiescent SMBHs are found to be correlated 
with the physical properties of their host bulges, such as the lu- 
minosity L (Kormendy & Richstone 1995; Marconi & Hunt 2003), 
the mass A/b (Magorrian et al. 1998; Marconi & Hunt 2003; Hiiring 
& Rix 2004), the stellar velocity dispersion a, (Ferrarese & Mer- 
ritt 2000; Gebhardt et al. 2000a; Tremaine et al. 2002; Ferrarese 
& Ford 2005), the galaxy light concentration Ckc (Graham et al. 
2001), the Sersic (1968) index n of surface brightness profile (Gra- 
ham & Driver 2007), and the iimer core radius r^, (Lauer et al. 
2007a). Among these, the Afbh-o"* relation of the form 

log(Afbh/M0) = Q + /31og(a/200kms-') (1) 

is one of the most tight with a intrinsic scatter of <0.3 dex (Novak 
et al. 2006). The validity of this relation have been explored for in- 
termediate mass black holes (IMBHs) in dwarf galaxies (Barth et 
al. 2005), the local active galaxies (McLure & Dunlop 2002; Green 
& Ho 2006), and Seyfert galaxies at 2 = 0.36 (Woo et al. 2006; 
Treu et al. 2007). In fact, this relation was predictd by Silk & Rees 
(1998), Fabian (1999) and Blandford (1999). As a crucial test of 
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theoretical models of formation and coevolution of the SMBHs and 
their host galaxies, many efforts have been made to explain the re- 
lation since its recognition (e.g. King 2003 and reference therein). 
The normalization, intrinsic scatter, and redshift-dependent evolu- 
tion of the relation are important for study of SMBHs demography 
(e.g. Yu & Tremaine 2002; Shankar et al. 2004). 

The early published estimation of the slope j3 are diverse from 
the highest value 4.8 ± 0.5 (Ferrarese & Merritt 2000, hereafter 
FMOO) to the lowest value 3.75 ± 0.3 (Gebhardt et al. 2000a, here- 
after GOO). Tremaine et al. (2002, hereafter T02) explored the rea- 
sons for the discrepancy between the results of the two groups, and 
show the systematic differences in the velocity dispersion measure- 
ments contribute the main diversity in the slope. The GOO and T02 
define "effective stellar velocity dispersion" cto as the luminosity- 
weighted rms dispersion within a slit aperture of length 2_Re, while 
FMOO use "central stellar velocity dispersion" CTc normalized to an 
aperture of radius Rc/8, where Rc is the effective radius of the 
galaxy bulges. T02 presented a best-fit result of the Afbh-o"o rela- 
tion with Q = 8.13 ± 0.06 and /3 = 4.02 ± 0.32, which is exten- 
sively used in the subsequent studies. The most current version of 
the A/bh-o-c relation are a = 8.22 ± 0.06 and jS = 4.86 ± 0.43 
(Ferrarese & Ford 2005, hereafter FF05). Although both of the re- 
lations are tight, their slopes still differ by 1.6 standard deviations. 

Besides the uncertainty in slope and normalization, the depen- 
dence of the Afbh-f* relation on the type of SMBHs host is also 
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worth consideration. A variety of observational and theoretical re- 
sults show that there are two kinds of "bulges" in disk galaxieQ 
(e.g. Kormendy & Kennicutt 2004, hereafter KK04). The classical 
bulges dwelling in the center of disk galaxies are like little ellip- 
tical galaxies, which are thought to be the products of violent re- 
laxation during major mergers. On the other hand, pseudobulges in 
disk galaxies are physically unrelated to ellipticals, their structure 
and kinematics resemble that of disks. They are believed to have 
formed via distinct scenario than that of classical bulges. There is 
no guarantee that SMBHs in pseudobulges should follow the same 
Mbh-o"* relation in ellipticals. In contrast, several evidences sug- 
gest they probably obey a different relation, (i) The rotationally 
flattened pseudobulges have smaller velocity dispersion than pre- 
dicted by the Faber- Jackson a,-L correlation (e.g. KK04). If the 
pseudobulge follow the same Alhh-L relation, they will naturally 
deviate from the same Mbh-fT", relation, (ii) The pseudobulges are 
dominantly shaped by the secular evolution, and probably never 
undergo a major merger, which is postulated as the key mechanism 
to explain the Mbh-o"* relation in the ellipticals. The gas inflow fu- 
eling SMBHs induced by the secular evolution is not as effective 
as that in major merger events, (iii) The slope of the Mhh-o', rela- 
tion in the local active galactic nuclei (AGNs) flattens somewhat at 
the low mass end (Greene & Ho 2006). Because most of the local 
low mass AGNs have late-type hosts (Kauffmann et al. 2003), we 
speculate some (or even most) of these low mass SMBHs may live 
in the pseudobulges, and disagree with the relation for high mass 
AGNs with (mostly) early-type hosts. 

Kormendy & Gebhardt (2001) checked the A/bh-o"* relation 
for five pseudobulges, and found no obvious inconsistency with the 
relation for the ellipticals. However, one of these galaxies (NGC 
4258) was in fact a classical bulge, incorrectly identified as a pseu- 
dobulge (J. Kormendy, private communication), and is consistent 
with the Mbh-o", relation for the elliptical galaxies, while the A/bh 
of the other four pseudobulges appears smaller than the predic- 
tion by the relation (cf. Fig. 5 of Kormendy & Gebhardt 2001). 
As pointed out by the authors, it is important to check their result 
with a larger sample. Therefore, in this paper we investigate the 
Afbh-o"* relation with the up-to-date sample of early-type bulges 
and pseudobuges. 

The paper is organized as follows. In section 2, we describe 
the sample used for analysis. In section 3, we determine and com- 
pare the A/bh-o"* relation in early-type bulges and pseudobulges. In 
section 4, we summarize and discuss our results. Throughout this 
paper, we use the the base 10 logarithms, and cosmological con- 
stants = 0.3, D.A = 0.7, h = 0.7. 



2 THE SAMPLE 

We select the sample with spatially well-resolved Afbh measure- 
ment from literature, and classify them as early-type bulges or 
pseudobulges. The selection criteria is 2rbh/fres >1, where rbh ~ 
GMbh/ff^ is the radius of the SMBHs influence sphere, rrcs is the 
instrumental spatial resolution. 

^ Athanassoula (2005) distinguished three types of bulges : classical bulges, 
boxy/peanut bugles and discy bulges. The boxy/peanut bugles are in fact a 
part of the bar seen edge-on; the discy bulges result from the inflow of 
gas and star formation in the central region of disk galaxies. Terminolog- 
ically, the "pseudobulges" in this paper include boxy/peanut bugles and 
discy bulges, as both of them are due to secular evolution and unrelated 
with the major mergers. 
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Figure 1. Comparison of the (Tc and ctc of the sample galaxies. The solid 
line denotes the model of T02 (eq. [2]). The dashed line denotes Cc = ctc- 

The properties of the bulges and SMBHs are listed in Table 1 
and Table 2. Our sample (39 early-type bulges -I- 9 pseudobulges) is 
larger than T02 (31 galaxies) and FF05 (30 galaxies), with updated 
mass and/or velocity dispersion measurement. Among the 48 ob- 
jects, 28 are weighted by stellar dynamics, 16 by gas kinematics (3 
have both stellar and gas kinematics measurements), 6 by masers, 
1 by stellar proper motions; 21 are inactive, 16 are LINERs, 2 are 
Seyfert 1, and 9 are Seyfert 2. 

Most distances of galaxies in the sample have been measured 
with the surface brightness fluctuation (SBF) method (Tonry et al. 
2001). The distances to the galaxies in Virgo cluster (NGC 4473, 
NGC 4486, NGC 4486A, NGC 4552, NGC 4564, NGC 4649) are 
based on the new SBF measurement by Mei et al. (2007). For those 
without an SBF distance we use recession velocities corrected for 
Virgocentric infall from HyperLeda databas^ except some indi- 
vidual galaxies explained in section 2.2. The black hole mass is 
linearly modified according to the distance that is different from 
the value used in the literature describing the mass measurement. 
The black hole mass errors are converted to Icr values if they are 
given at different confidence level (C.L.) in the literature. 

Most of the velocity dispersions CTc of E/SO galaxies are cal- 
culated using the formula a, = crco(8i?ap/i?c)° "'* (J0r gensen, 
Franx, & Kjaergaard 1995), where aco are provided by the Hy- 
perleda database, 7?ap ~ 0.85 kpc, Rc is the effective radius of 
bulges surface brightness profile, listed in the column (II) of Ta- 
ble 1 and 2. For Sa-Sd galaxies, observed a, profiles is not cen- 
tral peaked but nearly constant in the central 10" region, the varia- 
tion of (T, measured with different apertures is small (Pizzella et al. 
2004; Baes et al. 2003), so we take cr* value measured with an small 
aperture without further corrections. We adopt 5% error of (Jc and 
(Jc for E/SO galaxies and 10% error for Sa-Sd galaxies, except for 
some individual galaxies explained in section 2.2. For those with- 
out available Rc measurement, we simply use Rc = i?ap = 0.85 

^ http://leda.univ-lyonl.fr/ 
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Table 1. The sample of black holes in classical bulges and elliptical galaxies. 
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Notes: (1) name of the galaxy. (2) Hubble type of the galaxy. (3) the activity of the galactic nuclei (n = inactive; Syl = Seyfert 1; Sy2 = Seyfert 2; L 
= LINER). (4) distance to the galaxy. (5) mass of the central black hole. (6) method (s = stellar kinematics; g = gas kinematics; m = H2O masers; p = 
stellar proper motions) and reference for the black hole mass measurement. (7-9) effective and central stellar velocity dispersion of the bulge, and their 
references. (10-12) Sersic index and effective radius of the bulge surface brightness profile, and their references. 

" The value is taken from HyperLeda database, corrected with the aperture of Ra, the error is taken as 5% for the E/SO galaxies and 10% for the Sa-Sd 
galaxies. 

^ The value is estimated by the author, see details in section 2.2. 

The "core" elliptical galaxies (Lauer et al. 2007a). 
References: (1) Balcells et al. 2007; (2) Earth et al. 2001; (3) Batcheldor et al. 2005; (4) Bender et al. 2005; (5) Capetti et al. 2005; (6) Capetti & 
Balmaverde 2007; (7) Cappellari et al. 2002; (8) Cappellari et al. 2006; Cappellari et al. 2007a; (9) Cappellari et al. 2007b; (10) Copin et al. 2004; (11) 
Bower et al. 2001; (12) De Francesco et al. 2006; (13) De Francesco et al. 2008; (14) Devereux et al. 2003; (15) Kormendy et al. 1996; (16) Erwin et al. 
2002; (17) Faber et al. 1989; (18) Ferrarese et al. 1996; (19) Ferrarese & Ford 1999; (20) Gadotti 2008; (21) Gebhardt et al. 2000b; Shapiro et al. 2006; 
(22) Gebhardt et al. 2003; (23) Graham & Driver 2007; (24) Harms et al. 1994; Macchetto et al. 1997; (25) Neumayer et al. 2007; (26) Hermstein et 
al. 1999; (27) Hicks & Malkan 2008; (28) Houghton et al. 2006; (29) Kormendy & Bender 1999; (30) Laurikainen et al. 2004; (31) Laurikainen et al. 
2005; (31) Bower et al. 2001; (32) Marconi & Hunt 2003; (33) Nelson & Whittle 1995; (34) Nelson et al. 2004; (35) Nowak et al. 2007; (36) Onken et 
al. 2007; (37) Prugniel & Heraudeau 1998; (38) Richstone et al. 2004; (39) Saglia et al. 2000; (40) Sarzi et al. 2001; (41) Scarlata et al. 2004; (42) Silge 
et al. 2005; (43) Tremaine et al. 2002; (44) van der Marel & van den Bosch 1998; (45) Verolme et al. 2002; (46) Wegner et al. 2003; (47) Vega Beltran 
et al. 2001. 
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Table 2. The sample of black holes in pseudobulges. 
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Notes: Same as Table 1. 



^ The value is taken from HyperLeda database, corrected with the aperture of Rc (if available), the error is taken as 5% for the E/SO galaxies and 10% 
for the Sa-Sd galaxies. 

The mass uncertainty is not given in the literature, we use A log Mi)ii=0.3 as a rough estimation. 
References: (1) Bertola et al. 1995; (2) Davies et al. 2006; Hicks & Malkan 2007; (3) Drory & Fisher 2007; (4) Erwin et al. 2002; (5) Gadotti 2008; (6) 
Garcia-Rissmann et al. 2005; (7) Gebhardt et al. 2003; (8) Ghez et al. 2005; (9) Graham & Driver 2007; (10) Greenhill et al. 2003; (11) Kondratko et al. 
2008; (12) Lodato & Berlin 2003; Hure 2002; (13) Nelson & Whittle 1995; (14) Sarzi et al. 2001; (15) Terlevich et al. 1990; (16) Tilak et al. 2008; (17) 
Tremaine et al. 2002; (18) Yamauchi et al. 2004; Kondradco, Greenhill, & Moran 2005; (19) Fathi & Peletier 2003; (20) Shaw et al. 1993; (21) Maiohno 
et al. 1998. 



kpc, which is supposed to be an acceptable approximation to the 
true value. For those without Oc in literature, we take the value of 
o"c, since they are very close on average. T02 have modeled the 
(Jc-cTc relation for 39 early-type galaxies, and found (cf. eq. [16] in 
T02): 

log(^) = (-0.008±0.004)-(0.035±0.017)log(^^^^^^). 

(2) 

Comparison of the Oc and CTc of our sample galaxies is shown in 
Figure 1 . Obviously, the difference between the CTo and CTc are much 
smaller than their measurement errors. 

2.1 Identification of pseudobulges 

A comprehensive review on the pseudobulge properties and meth- 
ods to recognize pseudobulges are tabulated in Kormendy & Ken- 
nicutt (2004). Sersic index n < 2 is a simple but not safe criteria 
for pseudobulges (Fisher 2006), which is consistent with the obser- 
vations of our sample. Morphology and kinetic property are more 
reliable indicators than Sersic index. In general, if the "bulge" con- 
tains a nuclear bar, nuclear spirals, and/or nuclear ring, the bulge is 
identified as a pseudobulge. In (Vm/f7"b)-e diagram, pseudobulges 
locate above the oblate rotator line (e.g. Binney 1978), where Vm is 
the maximum line-of-sight rotational velocity of the bulge, (Tb is the 
rms stellar velocity dispersion of the bulge, and e is the character- 
istic ellipticity in the region interior to the radius of V^i- Bulges of 
different types (especially the boxy/peanut and discy bulges) often 
coexist in the same galaxy (Athanassoula 2005), the pseudobulges 
are found in all types (SO-Sm) of disk galaxies (Laurikainen et al. 
2007). The bulge properties of the disk galaxies in our sample are 
summarized in Table 3. 

NGC 1068 has nuclear spirals, it is identified as pseudobulges 
using HST optical images by Drory & Fisher (2007). 

NGC 2787 has a spectacular tilted nuclear dust ring of ra- 
dius ~1."5-10" (Sarzi et al. 2001; Erwin et al. 2003), a inner disk 
(Erwin 2003), and a bulge with disk-like kinetic property (Erwin 
2007). NGC 3384 has a inner disk (Erwin 2004) and a bulge with 
disk-like kinetic property (de Zeeuw et al. 2002). NGC 2787 and 



NGC 3384 are examples of the galaxies contain composite struc- 
tures consisting of both pseudobulges and small inner classical 
bulges (Erwin 2007). 

NGC 3079 has a peanut-shape bulge (Liitticke, Dettmar, & 
Pohlen 2000), which shows the cylindrical rotation (Shaw, Wilkin- 
son & Carter 1993), a character of disk-like pseudobulge. 

NGC 3227 has a asymmetric molecular nuclear ring with a 
diameter of ~ 3" (Schinnerer, Eckart, & Tacconi 2000), the Sersic 
index (n = 1.1) of the bulge is also disk-like (Gadotti 2008). 

NGC 3393 has nuclear spirals (Martini et al. 2003) and a pos- 
sible nuclear bar (Jungwiert et al. 1997), which identify the pseu- 
dobulge. 

The inner region of Circinus galaxy is bluer than typical bulge, 
that indicate young stellar populations (Wilson et al. 2000). The 
existence of a nuclear ring in CO emission (Curran et al. 1998) 
and near-UV light (Mufioz Marm et al. 2007) and a nuclear gas bar 
(Maiolino et al. 2000) identify the pseudobulge. 

IC 2560 has a boxy bulge (Lutticke, Dettmar, & Pohlen 2000). 
The red-shifted features in the H2O maser emission are slightly 
stronger than the blue-shifted, which is a typical characteristic seen 
in most megamasers and may be explained by a spiral structure of 
the maser-emitting nuclear disk (Ishihara et al. 2001). 

Milky way has a bar-related boxy bulge (e.g. Lopez- 
Corredoira et al. 2005), both the rotation (e.g. Minniti et al. 1992) 
and flattening (e.g. Sharpies et al. 1990) of the bulge suggest a pseu- 
dobulge. 

Most of the pseudobulges in our sample dwell in the barred 
disk galaxies (8/9), and most of our sample disk galaxies contain 
classical bulge are unbarred (11/14). Large scale bars are obvi- 
ously predominant component in disk galaxies with pseudobugles, 
and serve as a strong indicator for the pseudobugles. In the end, 
we caution that the identification of pseudobulges is still somewhat 
uncertain, some galaxies need further discrimination. 



2.2 Comments on individual objects 

Early-type bulges: 
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Table 3. Bulge properties. 



Galaxy 


Type 




n 




e 


Structures 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


N224 


Sb 


78 


2.2 






EB 


N1023 


SBO 


11 


2.0 


0.35 


0.33 


EB 


N3031 


Sab 


62 


2.0 








N3115 


SO 


82 


4.4 






BB 


N3245 


SO 


67 


4.0 






BB 


N3414 


SO 


11 


3.8 


0.09 


0.21 




N3998 


SO 


42 


4.1 








N4151 


Sa 


21 


3.0 








N4258 


SBbc 


71 


3.0 








N4459 


SO 


47 


3.5 








N4596 


SBO 


37 


1.4 


0.30 


0.15 




N5128 


SO 


48 


4.0 








N5252 


SO 


67 


6.0 








N7457 


SO 


76 


2.0 


0.62 


0.44 


BB 


N1068 


SBb 


33 


1.7 






NS 


N2787 


SBO 


66 


1.0 


0.91 


0.35 


NR, ID 


N3079 


SBcd 


86 








PB, NR 


N3227 


SBa 


69 


1.1 






NR 


N3384 


SBO 


59 


1.0 


0.77 


0.29 


ID 


N3393 


SBa 


31 








NS, NB? 


Circinus 


Sb 


63 








NR, NB 


IC 2560 


SBb 


66 








BB,NS? 


MW 


SBbc 


90 


1.3 






BB, ID 



Notes: (1) name of the galaxy. (2) Hubble type. (3) inclination of the disk 
are from Hyperleda or Ho 2007, except for N3384 (Neistein et al. 1999). 
(4) Sersic index of the bulge. (5) Vm / o"b ratio within Re - (d) ellipticity of 
the bulge. (7) EB = elliptical bulge, BB = boxy bulge, PB = peanut bulge 
(Liltticke. Dettmar, & Pohlen 2000); NR = nuclear ring; NB = nuclear bar; 
NS = nuclear spirals; ID = inner disk. 

NGC 224 (M31), NGC 3031 (M81) and NGC 3115 in our 
sample are typical classical bulges listed in KK04. 

NGC 221 (M32). — It is a small, low-luminosity elliptical, has 
a main body of Sersic profile with index n < 4 and extra light 
near the center above the inward extrapolation of the outer Sersic 
function (Kormendy 1999). 

NGC 1399.— The black hole mass (Mbh = 12l| x 10* Mq) 
is taken from Houghton et al. (2006), which is consistent with the 
resuk of Gebhardt et al. (2007), A/bh = (5.1 ± 0.7) x 10* Mq. 

NGC 3 1 1 5 and NGC 3245 have boxy bulges, while their bulge 
Sersic index is n ~ 4. We suggest central classical bulges and outer 
boxy bulges coexist in these galaxies. 

NGC 3377. — We estimate the Icr black hole mass error ac- 
cording to the Fig. 12 of Copin et al. (2004). 

NGC 3379. — The black hole mass is measure by stellar kine- 
matics (Afbh = l-OlJ:" X 10* Mq; Gebhardt et al. 2000b) and gas 
kinematics (Afbh = ^-^tl'l x 10* Mq; Shapiro et al. 2006), we 
take the average value of Mbh = (1.2 ± 0.6) x 10* Mq. 

NGC 3998. — The black hole mass is linearly modified with 
distance taken from Tonry et al. (2001), the la error is estimated 
from Fig. 12 of De Francesco et al. (2006). 

NGC 4151.— The black hole mass M^h = Stl x lO'' Mq is a 
mean value of the measurement using stellar kinematics (~4x 10^ 
Mq, Onken et al. 2007) and gas kinematics (3.2j;^;^ x lO'^ Mq 
[corrected for distance], Hicks & Malkan 2008), which is consis- 
tent with the RM measurement (4.57j;o 47 x 10^ Mq, Bentz et al. 
2006). The stellar velocity dispersion is taken from Nelson et al. 
(2004) measured with an aperture of 6". 5. 



NGC 4486A. — We estimate the la black hole mass error and 
velocity dispersion according to the Fig. 6 and Fig. 4 of Nowak et 
al. (2007). 

NGC 4596 has a relatively exponential bulge with Serisic in- 
dex n = 1.4 (Laurikainan et al. 2005), but the bulge kinematics 
(Bettoni & Galletta 1997) is ellipsoidal. 

NGC 4649. — We estimate the ae using the stellar velocity dis- 
persion profile measured by Bridges et al. (2006). 

NGC 5128.— The distance (D = 3.5 Mpc) and black hole 
mass (A/bh = 4.5to t x 10^ Mq) is taken from Neumayer et al. 
(2007), we estimate the la black hole mass error according to their 
Fig. 12. This black hole mass is consistent with the result of Haring- 
Neumayer et al. (2006), M^h = e.lt^J j x 10^ Mq. 

NGC 5252. — The stellar velocity dispersion is taken from 
Nelson & Whittle (1995) measured with an aperture of l".5 x 2". 3. 

NGC 7457. — This galaxy has a young box-shape bulge 
(Lutticke, Dettmar, & Pohlen 2000; Kormendy & Kennicutt 2004). 
Kormendy (1993) suggests it as a pseudobulge according to its 
small a*. But the kinetic property (Cappellari et al. 2007a) iden- 
tifies it as a classical bulge. 

The black hole masses of NGC 2974, NGC 3414, NGC 4552, 
NGC 4621, NGC 5813 and NGC 5846 are taken from Cappellari et 
al. (2007b) with a formal error of 10% (la C.L.), which is approx- 
imately converted from the 30% error at 3a" C.L. 

Pseudobulges: 

NGC 1068. — A prominent bar is present in the 2MASS in- 
frared image (Menendez-Delmestre et al. 2007), while the NED 
(RC3) type is Sb. The stellar velocity dispersion is taken from Nel- 
son & Whittle (1995) measured with an aperture of 1".5 x 2". 3. 

NGC 2787.— T02 gives a^ = 140 km s"\ while in the liter- 
ature much higher value are quoted, e.g., 210 km s^^ (Dalle Ore et 
al. 1991), 205 km s"^ (Neistein et al. 1999), 202 km s~^ (Barth et 
al. 2002), and 257 km s"^ (Erwin et al. 2003). As pointed by AUer 
et al. (2007), this discrepancy may stem from the central peak in 
the velocity dispersion profile. We estimate a^ and ac from Fig. 1 
ofBertolaetal. (1995). 

NGC 3079.— The accretion disk of the SMBH is probably 
thick, flared, and clumpy, in contrast to the compact, thin, warped, 
differentially rotating disk in the archetypal maser galaxy NGC 
4258 (Kondratko, Greenhill, & Moran 2005). The mass of SMBH 
thus may be underestimated. The stellar velocity dispersion is taken 
as the central value from Shaw et al. (1993). 

NGC 3227.— The mass 2.0t\i x Mq is a mean value 
of the measurement using stellar kinematics ([0.7-2.0] x 10^ Mq, 
Davies et al. 2006) and gas kinematics (2.31q 5 x 10^ Mq [cor- 
rected for distance]. Hicks & Malkan 2008), which is consistent 
with the RM measurement ([4.2 ±2.1] x 10^ Mq, Peterson et al. 
2004). The stellar velocity dispersion is estimated by the author 
from the Fig. 6, Fig. 1 1 and Fig. 16 of Davies et al. (2006). 

NGC 3393. — The stellar velocity dispersion is taken from Ter- 
levich et al. (1990) measured with an aperture of 2".l. 

Circinus. — The stellar velocity dispersion is taken as mean 
value of central 10" from Fig. 6 of Maiolino et al. (1998). The dis- 
tance is taken from Karachentsev et al. (2007), which is different 
from the value in Afbh measurement (Greenhill et al. 2003). 

IC 2560. — The stellar velocity dispersion is measured with an 
aperture of 2" x 2". 46 (Garcia-Rissmann et al. 2005). 

We do not include the following galaxies in our sample: 

NGC 1300 and NGC 2748.— NGC 1300 has spiral-like dust 
lanes down to the center, wrapping around the nucleus (Scarlata 
et al. 2004), and connecting to the leading edge of a strong kpc- 
scale bar (Martini et al. 2003), which suggests gas inflow in both 
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large and small scale. The bulge of NGC 1300 is disk-like with 
Sersic index n — 1.3 (Laurikainen et al. 2004), and identified as a 
pseudobulge by Kormendy & Fisher (2005). The nuclei of N2748 
is not well defined in HST observations due to dust lane obscura- 
tion (Scarlata et al. 2004). The masses of SMBHs in NGC 1300 
(6.6t^i X 10^ Mq) and NGC 2748 (4.4i|J x lO'^ M©) are mea- 
sured by nuclear gas kinematics (Atkinson et al. 2005). The authors 
neglected the effect of extinction by dust in the analysis of the 
gaseous rotation curves. The extinction in the H band which was 
used to measure the black hole masses yield an underestimation in 
luminosity of ~ 20% for NGC 1300 and ~ 58% for NGC 2748 in 
their central 50 pc regions. Furthermore, the extinction is probably 
uneven due to the complex dust structures. Therefore, the A/bh in 
two galaxies are not well estimated. Their results are much higher 
than the prediction of Mbh-o"* or Mhh-L relations (7 x 10^ Mq or 
2 X 10'^ Mq for NGC 1300; 6 x lO'^ Mq or 7 x 10*' Mq for NGC 
2748). In fact, the HST can only barely resolved the sphere of influ- 
ence of SMBHs with the mass predicted by the A/bh-o"* or Alhh-L 
relations (rbh = 0".04 ~ 0".l for NGC 1300; rbh = 0".03 for 
NGC 2748). 

NGC 2778 and NGC 4303.— //5rcannot resolve the influence 
radius of the SMBHs in NGC 2778 (2rbh/rres = 0.8, Gebhardt et 
al. 2003) and in NGC 4303 (2rbh/rros = 0.9, Pastorini et al. 2007). 
The zero black hole mass model for NGC 2778 is ruled out at only 
the 95% C.L. (Gebhardt et al. 2003). 

NGC 4342. — The mass of the central SMBH measured by 
Cretton & van den Bosch (1999) make this object a outlier for the 
Mhh-L relation (Marconi & Hunt 2003) and the A/bh-A/b relation 
(Haring & Rix 2004). According to the analysis of Valluri, Merritt 
& Emsellem (2004), no best-fit value of Afbh can be found for this 
object, the published value 3.0'^\'q x 10* M© is probably overes- 
timated. 

NGC 4350 (Pignatelli et al. 2001), NGC 4486B (Kormendy 
et al. 1997), and NGC 4594 (Kormendy et al. 1996).~Their mass 
measurement is uncertain due to lack of three-integral stellar dy- 
namical models. 

NGC 4374. — As pointed by T02, the published measurement 
of A/bh in NGC 4374 (M84) is differ from each other beyond the 
stated errors. 

NGC 7582.— The influence radius of the SMBH is not well 
resolved, 2rbh/rrcs = 0.7 (Wold et al. 2006). 

The more reliable and accurate mass measurement of these 
SMBHs in the future will contribute to the poorly sampled low- 
mass end of the A/bh-o"* relation, especially for the pseudobulges. 



3 THE ESTIMATION OF A/bh-ct, RELATIONS 
3.1 The linear fitting algoritlim 

Different authors use different algorithms to estimate the log-linear 
Afbh-o"* relation eq. ( 1 ). GOO and T02 fit the relation by minimizing 
"X^ estimator" (Press et al. 1992) 

2 ^ (y. - g - f3x,) 

where x — log(cr/200 km s""*^), y = log(Afbh/M0), e^i and tyi 
are errors of x and y measurements. In order to account for the 
intrinsic scatter in the Afbh-o"*, T02 replace eyi by (e^i + t^Y^^ , 
where eo is a constant denotes the intrinsic scatter of j/ at a fixed 
X. The value of eo is set so that the best-fit reduced is equal to 
its expectation value of unity, T02 find eo = 0.27 for their sample. 



The 1(7, 2(7, 3a uncertainties in a and (5 are given by the upper and 
lower limits of a and f) for which ~ Xmin — 1-01, 6.75, 16.68. 

FF05 use the "AB estimator" described by Akritas & Bershady 
(1996) to determine the best-fit a, [3 and their uncertainties. Both 
and AB estimators have the advantages of accounting for mea- 
surement uncertainties in both variables and are asymptotically nor- 
mal and consistent. In our analysis, the AB estimator gives the sim- 
ilar estimate for the slope /? (cf. Figure 2, Figure 3, and Table 4). A 
review of the comparison of the two estimator can be found in T02. 

In the subsequent analysis, we use both logarithms to estimate 
a and /?, and denote them as a^2 , (i^i and OfAB, /3ab, respectively. 



3.2 Tlie slope and normalization of tlie relations 

The fitting results of the Afbh-o"c relation and the Afbh-o"c relation 
are shown in Figure 2. The best-fit value and uncertainties of q, 
P and eo are summarized in Table 4. We notice the Mhh-o'c and 
Afbh-o"c relations are roughly the same, and the results derived by 
two fitting algorithms are very close, both within Icr uncertainty. 

The slope /3 for the early-type bulges derived by two fitting al- 
gorithms and two dispersion definitions, is in the range of 3.9~4.1 
with uncertainties of <0.28. It is consistent with the result of T02 
(/3 = 4.02 ±0.32). We take the mean fitting results by estimator 
as the "uniform" A/bh-o"* relation for the early-type bulges: 

a = 8.28 ±0.05, /3 = 4.06 ± 0.28. (4) 

The case for the pseudobulges are similar, but with the different 
value. Their best-fit of A/bh-o", by x^ estimator is 

Q = 7.50 ±0.18, /3 = 4.5±1.3. (5) 

The intrinsic scatter for the early-type bulges and pseudobulges are 
0.27 and 0.25 respectively. The Afbh measurement and the predic- 
tion by the best-fit relation are compared in the middle panels of 
Figure 2. The best-fits of T02 and FF05 are plotted in the same 
figures for comparison. Most of the pseudobulges lie below the 
prediction of relation for the early-type bulges and of all the three 
previous formula. 

The discrepancy between the A/bh-o", relations for the early- 
type bulges and pseudobulges is explicit in the bottom panels 
of Figure 2. They are distinct over 3(7 significance level. The 
Kolmogorov-Smimov test gives a very low probability (P < 2 x 
10~*) that the Afbh-cr* relation for pseudobulge agrees with the 
relation for the early-type bulges. 

The relation for the pseudobulges has a similar slope but 
0.78 ± 0.18 dex smaller normalization than the relation for the 
early-type bulges. On average, the SMBHs in early-type bulges are 
as 6 times massive as ones in pseudobulges with the same velocity 
dispersion. T02's normalization (a = 8.13 ± 0.06) is just between 
our estimates of the early-type bulges and pseudobulges, may be a 
result of fitting for a mixed sample of two types of bulges. 

We also fit the A/bh -(7* relation for a combined sample of all 
the objects, obtain a steeper slope than that of the both samples. 
It is easy to be understood that, the A/bh-<7, relation for pseudob- 
ulges is roughly parallel but lower than that for early-type bulges, 
the high mass end dominant early-type bulges and the low mass 
end dominant pseudobulges will thus determine a steeper relation 
for the combined sample. It can partially explain the steep slope 
derived by FF05 (cf. the right panels of Figure 2). 
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Figure 2. Left-top: The Af^h-o'c relations for early-type bulges (filled markers, thick lines) and pseudobulges (open markers, thin lines). Mass measurement 
based on stellar kinematics are denoted by squares, on gas kinematics by circles, on masers by triangles, and on stellar proper motions by diamond. The solid 
lines are fitting results using estimator, dashed lines using AB estimator. Left-middle: Residuals between the measured A/bh and the prediction by best-fit 
relation using estimator for early-type bulges. The markers are the same as the top panel. The best-fit relation of T02 is shown for comparison. Left-bottom: 
The contours are Icr, 2ct and 3(t limits of x^ fitting of a and /3. The central circles denote best-fit value. The dashed squares and their central cross markers 
denote the best-fit value and uncertainties using AB estimator. The thick and thin lines are results of early-type bulges and pseudobulges respectively. Right: 
Similar as the left panel, but for the Afbh-Cc relations, the best-fit of T02 is substituted by one of FF05. 
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Table 4. The fitting results of the Afbh"""* relation. 



















A^bh-O'c 






Host type 






£0 


"AB 


/3ab 






eo 


"AB 


/3ab 


early-type bulges 


8.30±0.05 


4.10±0.28 


0.26 


8.31±0.05 


4.07±0.24 


8.26±0.05 


4.01 ±0.28 


0.27 


8.27±0.05 


3.97±0.22 


pseudobulges 


7.50±0.18 


4.46±1.30 


0.25 


7.51±0.10 


4.25±0.41 


7.51±0.17 


4.62±1.32 


0.24 


7.50±0.10 


4.20±0.40 


all galaxies 


8.18±0.06 


4.57±0.37 


0.42 


8.18±0.06 


4.50±0.30 


8.14±0.06 


4.52±0.36 


0.39 


8.14±0.06 


4.59±0.32 


core ellipticals 


8.28±0.11 


4.84±0.78 


0.24 


8.27±0.16 


5.62±1.47 


8.17±0.11 


5.33±0.81 


0.21 


8.15±0.21 


6.25±1.79 


normal ellipticals 


8.28±0.07 


3.97±0.38 


0.29 


8.28±0.07 


4.00±0.36 


8.22±0.06 


3.82±0.36 


0.29 


8.23±0.07 


3.81±0.31 




Figure 3. Left-top: The best-fit Mth-Cc relations for "core" galaxies (open circles, thin lines) and the other early-type galaxies (filled circles, thick lines). The 
solid lines and dashed lines denote fitting results using ™d AB estimator respectively. Left-bottom: The contours are Icr, 2rj and 3(t limits of fitting of 
a and /3 in eq. (1). The central circles denote best-fit value. The dashed squares and their central cross markers denote the best-fit value and uncertainties using 
AB estimator The thin and thick lines are results of the "core" elliptical galaxies and the normal elliptical galaxies. Right: Similar as the left panel, but for the 
A/bh-fc relations. 



3.3 The Mbh-cr, relation for the "core" elliptical galaxies 

Some elliptical galaxies are found to have cores in their central 
brightness profiles (Lauer et al. 1995). A core is evident as a radius 
at which the steep envelope of the galaxy breaks and becomes a 
shallow profile in logarithmic coordinates. The cores are believed 
formed by binary BHs in the dissipationless ("dry") mergers. These 
remnants of dry mergers will preserve the Mbh-i relation but devi- 
ate from the Mbh-o"* relation, because the Faber- Jackson relation- 



ship L oc for the "core" galaxies are much steeper than than 
the classic L oc cr* for the "normal" elliptical galaxies (Lauer et al. 
2007a). On the other hand, the SMBHs masses Mbh predicted by 
the i\fbh-cr« relation is systematically smaller the masses predicted 
by the Mbh-i relation for high-luminosity galaxies, such as bright- 
est cluster galaxies (BCGs). The BCGs have large cores, which are 
probably formed in the dry mergers. Hence, the Mbh-o"* relation 
will break at the high mass range due to the core galaxies made by 
dry mergers. 
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Inspired by these arguments, we attempt to compare the Mbh- 
(j, relations of core elliptical galaxies in our sample and of the other 
normal early type galaxies. We augment the 13 core elliptical galax- 
ies in Table 1 with the brightest group galaxy (BGG) Cyg A (Tad- 
hunter et al. 2003; Thornton et al. 1999), and the BCGs in galaxy 
cluster Abell 1836 and Abell 3565 (Dalla Bonta et al. 2007). The 
BGGs and BCGs are believed to have similar physical origin with 
core elliptical galaxies. Their central black hole mass measured by 
gas kinematics and central stellar velocity dispersiorPlare: Cyg A, 
Mbh = [2.5 ± 0.7] X lO'' Mq, (j,=290±70 km s~^ Abell 1836 
BCG, Mbh = 4.8tg j X 10^ M0, cr.=300±50 km s^^ Abell 3565 
BCG, Mbh = l-SlJ^i X 10^ Mq, cr.=200±50 km s"\ Figure 3 
depict the A/bh-c* relations for the core elliptical galaxies and the 
other normal elliptical galaxies. The best-fit results are presented in 
Table 4. Two relations are still tight (eo J; 0.3), consistent within 
la uncertainties, but the former is steeper (/3 ~ 5-6) than that later 
(P ~ 4). The mean values of the fitting results by estimator are 

a = 8.23 ±0.11, /3 = 5.09 ±0.80, (6) 

while the AB estimator gives a steeper slope (/3 ~6). 

Our result confirms the prediction of Boylan-Kolchin et al. 
(2006) that the i\/bh-c"« relation for the core elliptical galaxies is 
steeper than that for the normal galaxies. The Mbh-cr* relation de- 
rived from all the elliptical galaxies underestimates the mass of 
SMBHs in the center of core galaxies statistically. For the most 
massive galaxies (a, ~ 400 km s~^), the discrepancy can be 
as large as factor 2~3. The sub-sample of core galaxies over- 
lap the one of the normal elliptical galaxies in the mass range of 
10* ~ 10^'^ M0, therefore the intrinsic scatter of the correlation 
for the combined sample in the high mass range is larger than that 
in the low mass range. 



4 CONCLUSIONS AND DISCUSSION 

We examine whether early-type bulges and pseudobulges obey the 
same relation between black hole mass and stellar velocity disper- 
sion, using the largest published sample to date. The main conclu- 
sions of this paper are as follows: 

(1) On the relation for early-type bulges, we find a simi- 
lar slope (P — 4.06 ± 0.28) and higher normalization (a — 
8.28 ± 0.05) than the previous result of Tremaine et al. (2002). 
The intrinsic scatter of the relation is <0.27 dex. 

(2) The pseudobulges have their own tight Afbh-o", relation 
with /3 = 4.5 ± 1.3, a = 7.50 ± 0.18, which are different from the 
relation for early-type bulges over the 3a significance level. 

(3) The "core" elliptical galaxies at high black hole masses 
end follow a steeper A/bh-f* relation with slope /3 ~ 5-6. 

The mass density of local SMBHs calculated with the A/bh-o"« 
relation (e.g. Yu & Tremaine 2002; Wyithe & Loeb 2003; Shankar 
et al. 2004) is consistent with the cumulative mass density of optical 
or X-ray selected quasar remnants expected from energy arguments 
(e.g. Soltan 1982; Fabian & Iwasawa 1999; Yu & Tremaine 2002; 
Shankar et al. 2004; Marconi et al. 2004). However, the 0.15±0.08 
dex higher normalization of our best-fit relation for the early-type 
bulges than the value of T02 will enhance the predicted mass and 
mass density of local SMBHs by ~ 40%. In addition, the steeper 
slope of the relation for the core galaxies will raise the mass func- 
tion of the most massive black holes in the universe, while the lower 

^ These objects are excluded in section 2 for their relatively uncertain cr, . 
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Figure 4. Comparison of the Afbh'O'o relation for elliptical galaxies (stars), 
classical bulges (filled squares), and pseudobulges (open squares). The 
ban'ed disk galaxies are marked by circles. The thick and thin solid lines 
are the best fit results for the early-type bulges and the pseudobulges re- 
spectively. 

relation for the pseudobulges will reduce the mass function of the 
low mass SMBHs in the centers of disk galaxies. Such modifica- 
tion will also affect estimation of the fraction of obscured AGNs, 
which can be tested independently by the observations of cosmic 
hard X-ray background. We will elaborate these implications in the 
subsequent papers. 

We note unbiased estimate for parameters /3 and a requires 
larger unbiased sample (both in A/bh and a,), or one must deter- 
mine the selection bias carefully (e.g. Lauer et al. 2007b). However, 
the significant lower normalization (~ 0.8 dex) of the Mbh-o"* re- 
lation for pseudobulges seems not due to the effect of selection 
bias. At a given velocity dispersion, the SMBHs with massive Ai^bh 
are easier to be discovered and measured, because their influence 
spheres is larger. Our findings can be tested in larger unbiased sam- 
ples, especially for the pseudobulges. The preliminary results of the 
Afbh upper limit on three SMBHs in pseudobulgefj (NGC 3351, 
NGC 4143, NGC 4314) indicate they also he below the Ahh-a* 
relation for early-type bulges (Corsini et al. 2007). 

Another interesting and unexpected finding is the agreement 
between the slope of Afbh-o"G relation and of the A/bh-o"c rela- 
tion. T02 concluded that the earlier found discrepancy arises mostly 
from the different definitions of stellar velocity dispersion, but this 
difference in our sample is insignificant (cf. Figure 1). Our esti- 
mate of ac is much shallower than the previous results of FMOO 
and FF05. This "converged" slope will be important, and in some 
senses, convenient for the theoretical models and numerical simu- 
lations that try to reproduce the relation. 

We compared the A/bh-o"* relation for elliptical galaxies, clas- 
sical bulges, pseudobulges, and barred/barless disk galaxies. As 
shown in Figure 4, elliptical galaxies and classical bulges follow 
the same A/bh-o"» relation. In our sample, barred/barless distinc- 
tion is only a result of classical/pseudobulges dichotomy. We also 

They ai'e identified as pseudobulges by the nuclear rings for NGC 3351 
and NGC 43 14 (Shaw et al. 1995), and by the nuclear spirals for NGC 4143 
(Martini et al. 2003). 
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Figure 5. Comparison of the inclination distribution of tlie disk galaxies 
with classical bulges (open columns) and pseudobulges (grey columns). 



compared the inclination of the disk galaxies with classical bulges 
and pseudobulges (Figure 5). No obviously different distribution is 
found in our sample, the inclination seems not be responsible for 
the different Afbh-c"* relations. However, we need larger sample to 
consider the effect of bar and inclination on the Mbh-o"* relation. 

In stead of the assumed universal A/bh-o"» relation, the real 
correlation is composed of three sub-samples of SMBHs: hosted 
by pseudobulges, by normal elliptical galaxies or classical bulges, 
and by "core" elliptical galaxies. They have different slope and 
normalization, but similar intrinsic dispersion. The tight correla- 
tion between A/bh and cr, demonstrate that the bulge formation and 
SMBHs fueling are closely connected, despite the type of bulges. 
The strong feedback due to wet major mergers is the favored mech- 
anism responsible for the formation of the correlation for the early- 
type bulges, while the slight steeper correlation for the core galax- 
ies is build by the dry mergers of their early type progenitors. If 
the SMBHs in pseudobulges grow in the same period of their host 
bulge formation, their tight correlation may be the result of the 
secular evolution. The similar slope of the relations for early-type 
bulges and pseudobulges indicates they may experience the simi- 
lar AGN feedback processes (e.g. Silk & Rees 1998; King 2003), 
but the SMBHs growth in pseudobulges is relatively slower than 
that in the early-type bulges. We suggest it is due to the lower effi- 
ciency of AGN fueling by secular evolution. The other possibility 
is that the formation of these SMBHs is earlier than their host pseu- 
dobulges. The newborn disk galaxies may have the same A/bh-o"* 
relation with that of elliptical galaxies, but the secular evolution 
steadily enhanced the central velocity dispersion, and the subse- 
quent growth of SMBHs is insignificant. Therefore, the pseudob- 
ulges move rightward in the A/bh-cr, diagrain and deviate from the 
original relation. In conclusion, the SMBHs formation and evolu- 
tion scenario in the pseudobulges deserve further exploration, es- 
pecially for the early stage of disk galaxy formation. 
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